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ABSTRACT 

We use an invariant operator to study the quantum gravitational back- 
reaction to scalar perturbations during chaotic inflation. Our operator is the 
inverse covariant d'Alembertian expressed as a function of the local value of 
the inflaton. In the slow roll approximation this observable gives —1/(2H 2 ) 
for an arbitrary homogeneous and isotropic geometry, hence it is a good 
candidate for measuring the local expansion rate even when the spacetime is 
not perfectly homogeneous and isotropic. Corrections quadratic in the scalar 
creation and annihilation operators of the initial value surface are included 
using the slow-roll and long wavelength approximations. The result is that 
all terms which could produce a significant secular back-reaction cancel from 
the operator, before one even takes its expectation value. Although it is not 
relevant to the current study, we also develop a formalism for using stochastic 
samples to study back-reaction. 
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1 Introduction 



Cosmological perturbations generated during an inflationary era are almost 
certainly responsible for the density ripples out of which the presently ob- 
served stars, galaxies and clusters have grown. The basic inflationary mecha- 
nism that generated the seeds for galaxy formation and growth of large-scale 
structure is quantum mechanical particle production powered by the accel- 
erated expansion of the Universe. Now that the latest balloon experiments 
[[!]] have all but validated this picture, we are left with the duty of testing 
further the consequences of particle production in the early Universe. An 
outstanding question is whether quantum pair creation has any effect on the 
background in which it takes place. 

The simplest system in which the problem of back-reaction can be studied 
is scalar-driven inflation. Since dynamical scalars mix with the gravitational 
potentials they cause, gravitational back-reaction can in principle occur at 
one loop in scalar-driven inflation. Indeed, previous results by the present 
authors and collaborators 0, §, f|, || suggested that the one loop back-reaction 
of scalar perturbations could slow down the expansion rate of the Universe 
in the case of the so-called "chaotic" inflationary models. 

In RJ, spatial averaging and a fixed gauge were employed in order to 
compute the effective energy-momentum tensor for cosmological perturba- 
tions. The present authors studied the same physical problem |4|, ||, this 
time taking expectation values of the metric and using those expectation 
values to form invariants, and came to conclusions identical to those of @, §]. 
Nevertheless, these works have been criticized in two grounds. 

The first objection || concerns the use of expectation values of the gauge- 
fixed metric to build physical observables. The concern is that certain vari- 
ations in the gauge fixing condition could change the expectation value of 
the metric in ways which cannot be subsumed into a coordinate transforma- 
tion. Therefore, forming the expectation value of the metric into coordinate 
invariant quantities would not purge these quantities of gauge dependence. 
Back-reaction should then be studied with an operator which is itself an 
invariant, before taking the expectation value. 

The second objection ]/]] states that using expectation values invites a 
Schrodinger Cat paradox. This is because the process of superadiabatic am- 
plification leaves superhorizon modes in highly squeezed states whose behav- 
ior is essentially classical. That is, the various portions of the wave function 
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no longer interfere with one another. Quantum mechanics determines the 
random choice of where we are in the wave function, but the evolution of 
each portion is approximately classical. Under these conditions, averaging 
over the full wave function does not give a good representation of the physics. 
In the context of the Schrodinger Cat paradox Linde observes: "When you 
open the box a week later what you find is either a very hungry cat or else a 
smelly piece of meat, not the average of the two." He and his colleagues would 
prefer that back- reaction be studied stochastically ||, where each mode is 
assigned a random, (D-number value as it experiences horizon crossing and 
evolves classically thereafter. 

In the present study of the one-loop back-reaction effect we have tried to 
address these objections. To avoid potential problems from using the gauge 
fixed metric we have instead computed the functional inverse of the covariant 
d'Alembertian (see also [jj for a more detailed discussion): 




This operator — A = □ c 7 1 [g](i, x) — averages over the past lightcone, very 
much as astronomers do when compiling a Hubble diagram. In the slow roll 
approximation it gives —\H~ 2 for an arbitrary homogeneous and isotropic 
universe ||. The scalar A is therefore a reasonable candidate for measuring 
the local expansion rate when the universe is not precisely homogeneous and 
isotropic. By measuring A on 3-surfaces upon which the dynamical scalar is 
constant we promote it into a full invariant, A\ nv . 

At one loop, and to leading order in the standard infrared expansion, 
corrections to the operator A depend upon a variable we call This is 

the accumulated Newtonian potential from infrared modes. <&(x) grows with 
time, since more and more modes redshift into the infrared regime as inflation 
progresses, so dependence upon any positive power of would give a 

secular effect. However, when full invariance is enforced by evaluating A on 
3-surfaces of constant scalar, all terms involving only the undifferentiated 
Newtonian potential cancel. This result was anticipated by Unruh || who 
noted that the linearized mode functions become pure gauge in the long 
wavelength limit. The fact that more and more modes redshift to the infrared 
regime is not pure gauge, but Unruh's observation means that physical effects 
must involve derivatives of Such terms do contribute, but they are 
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small. We conclude that there is no significant deviation from the background 
expansion rate at this order, in marked contrast with the results of || [|, [|]. 

Section 2 describes the perturbative background in which we work. It also 
illustrates the slow roll expansion through which one can perform most oper- 
ations of temporal calculus in an arbitrary inflationary background. Section 
3 describes the formalism by which the field equations are solved perturba- 
tively. It concludes with the long wavelength approximation which effectively 
removes spatial dependence and produces a one dimensional problem that 
can be treated in the slow roll expansion. Section 4 applies the technology 
to determine the various field operators to quadratic order in the creation 
and annihilation operators on the initial value surface. The operator A is 
evaluated to the same order in Section 5. The distinction between expec- 
tation values and stochastic samples is irrelevant, at this order, in view of 
the cancellation of all potentially significant one loop corrections to the A- inv 
operator. However, the issue should re-emerge at higher loops where secu- 
lar back-reaction can derive from the coherent superposition of interactions 
within the observer's past lightcone. We have therefore developed a formal- 
ism for the perturbative study of stochastic effects, which is presented in 
Section 6. Our results are summarized and discussed in Section 7. 

2 The perturbative background 

The system under study is that of general relativity with a general, minimally 
coupled scalar: 

C = Y^g R ^ ~ \d^d v ^^~g - V{ V )^~g . (2) 

This section concerns the homogeneous and isotropic backgrounds go and 
ipo about which perturbation theory will be formulated. Three classes of 
identities turn out to be interesting for our purposes: 

1. Those which are exact and valid for any potential V(ip); 

2. Those which are valid in the slow roll approximation but still for any 
potential; and 

3. Those which are valid for the slow roll approximation with the potential 
V(<p) = |mV- 
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We shall develop them in this order, identifying the point at which each 
further specialization and approximation is made. 

Among the exact identities is the relation between co-moving and confor- 
mal coordinates: 

dsl = —dt 2 + a 2 dx ■ dx = a 2 (rj) j— drj 2 + dx ■ dx} , (3) 

where conformal time r] is defined in the usual way: 

dt = a{rj) dr] . (4) 

The Hubble "constant" is the logarithmic co-moving time derivative of the 
background scale factor: 

H^-= a -, (5) 

a a 2 

where a dot denotes differentiation with respect to (background) co-moving 
time and a prime stands for differentiation with respect to conformal time. 
Two of Einstein's equations are nontrivial in this background: 

3# 2 = \^{\^l + V{^)) , (6) 
-2H-3H 2 = ^ 2 {^-%)} , (7) 

where k 2 = 16nG is the loop counting parameter of perturbative quantum 
gravity. By adding the two Einstein equation one can solve for the time 
derivative of the scalar, which we assume to be negative: 



<Po = —yJ-H. (8) 

K 

Successful models of inflation require the following two conditions which 
define the slow roll approximation: 

\<p \ < if |£o| , (9) 
0o « V(<p ). (10) 

It follows that there are two small parameters. Although these are tradition- 
ally expressed as ratios of the potential and its derivatives the more useful 
quantities for our work are ratios of the Hubble constant and its derivatives: 

^«i , (11) 

H 2 -HH 
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For models of interest to us the rightmost of these parameters is negligible 
with respect to the leftmost one. 

The slow roll approximation gives useful expansions for simple calculus 
operations. For example, ratios of derivatives of the field are: 




(12) 



f = 2 ^(l + ^ + ...J. ,13) 
Successive partial integration also defines useful slow roll expansions: 



J dtH a a 
J dtH a 




(15) 



3 Einstein-scalar during inflation 

The purpose of this section is to describe our procedure for expressing the 
metric and scalar fields at any spacetime point as functionals of the uncon- 
strained initial value operators. These unconstrained initial value operators 
are the only true degrees of freedom in any quantum theory and expressing 
the dynamical variables as functionals of them is what it means to solve the 
equations of motion in the Heisenberg picture. As a simple example, if the 
dynamical variable is the position q(t) of a one-dimensional particle then the 
equation of motion, 

q(t)+u 2 q(t) = 0, (16) 

has as its "solution" , 

?(0 = Qo cos(c<jt) + — sin(cjt) . (17) 

In this case the unconstrained initial value operators are go and go. 

Obtaining such a solution perturbatively entails two expansions. In the 
first stage one writes the dynamical variables as background plus perturba- 
tions and expands the equations of motion in powers of the perturbed quan- 
tities. To continue with our particle example, suppose the full equation of 
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motion is = and that we are perturbing around some classical solu- 

tion X{t). If we define the perturbed position as q(t) then x(t) = X(t) + q(t). 
Since the background is a solution the perturbative equations of motion nec- 
essarily begin at linear order, 

E[X + q](t) = q(t) + uj 2 q(t) + AE[q](t) , (18) 

where AE[q](t) contains terms of order q 2 and higher. 

The second stage of perturbation theory consists of solving for the per- 
turbed quantities as a series in powers of the unconstrained initial value oper- 
ators. This is done by first solving the linearized equations so as to make the 
unconstrained variables agree with their full initial value data on the initial 
value surface. The linearized solution in our particle example would be, 



Sn„(t) = ^{ae-' + ate^} , (19) 

where we have chosen to organize the initial value data into the annihilation 
operator, a = muj / 2h(q + iq /uj), and its conjugate. 

One then "integrates" the full equations of motion using the retarded 
Green's functions of the linearized theory. The retarded Green's function of 
our particle example is, 

G Ict (t; t') = 9 JLlIl sin [^ _ t ')] , (20) 
uj 

and the integrated equation of motion is, 

q(t) = qun(t) - - f df sinK* - t')]AE[q}(t') . (21) 
UJ Jo 

The perturbative solution derives from successive substitution of this inte- 
grated form for the perturbed quantities on the righthand side of the equa- 
tion, 

q(t) = qUt) - - f dt' sm[u(t - t')]AE [q hn ] (if) + ... . (22) 
uj Jo 

It will be seen that the perturbative operator expansion which results 
bears a close relation to the diagrammatic expansion of quantum field the- 
oretic amplitudes. The place of external lines is taken by the linearized 
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solution; propagators are replaced by the corresponding retarded Green's 
functions; and, except for some factors of i, there is no change at all in the 
vertices. Since no loop diagrams appear there is no need of Faddeev-Popov 
ghosts. Owing to this correspondence many of the results we shall obtain for 
the operator expansion have already appeared in our previous perturbative 
expansion for expectation values J| . We shall retain most of the conventions 
used in that paper. 

In a gauge theory such as gravity there is the additional complication of 
constraints which require some of the dynamical variables to be nonlinear 
functionals of the unconstrained initial value operators even on the initial 
value surface. We shall deal with this by working in a gauge with a nonlocal 
field redefinition such that there are no nonlinear corrections on the initial 
value surface. Of course the nonlinear corrections required by the initial 
value constraints return when our field redefinition is inverted to give the 
perturbed metric and scalar fields. 

Another complication is that one cannot generally obtain explicit solu- 
tions to the linearized field equations in the time dependent backgrounds of 
chaotic inflation. We therefore explain how to construct plane wave solutions 
as series expansions around the ultraviolet and infrared limiting cases. Since 
the putative physical effect derives from modes which have redshifted into 
the infrared regime it is the long wavelength approximation — coupled with 
the slow roll expansion — which allows us to obtain quantitative results. 

Finally, it should be noted that we are not actually including all the 
physical degrees of freedom. Since gravitons cannot give a significant back- 
reaction at one loop we have suppressed the initial value operators associated 
with them, leaving only scalar degrees of freedom. Of course we still work 
out how the full metric depends upon this scalar initial value data. 

3.1 Perturbative field equations 

The field equations derived from the Einstein-scalar action ([2]) are most use- 
fully expressed in the form, 



F 




(23) 



a 




. 



(24) 
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Here G^ u = R^ u — -g^R is the Einstein tensor and T^ u is the scalar stress 
tensor whose covariant expression has the form, 

Tfu, = dptpdvip - g^ (^g^dpcpd^ + V(<p^j . (25) 

Our dynamical variables are the scalar perturbation <p and the conformally 
rescaled pseudo-graviton ip^: 

(p(r),x) = (p (rj) +0(?7, x) (26) 
g^iViX) = a 2 (r]) [rj^ + Kip^t], x)} , (27) 

where r]^ is the (spacelike) Lorentz metric which is used to raise and lower 
pseudo-graviton indices. We define the quantities AF and AE^ IU to include 
all terms of second and higher orders in ip^ and <p. One can read off the 
quadratic terms of AF and AE^ U from Tables 1 and 3 of [§J: multiply the 
terms of Table 1 by 1 /a and then drop the "external" ip^ v ; vary the terms of 
Table 3 with respect to (p or ip^ and then multiply by 1/a. 

The linearized equations are vastly simplified by imposing the gauge con- 
dition, 

F^ = a (^p^ u - \^,n - 2 ~^M> + ^o'Vo^ = , (28) 

where a comma denotes ordinary differentiation and ip = Tj^tp^. It is also 
useful to 3 + 1 decompose and to slightly rearrange the perturbed fields as 
follows, 

/ = a<p , z = aip 00 , Vi = aip 0i , h i5 = a (^ - c^Voo) • (29) 

The various differential operators of the linearized equations can also be given 
a simple common form, 

Vj = V 2 + Vj = V 2 - dl + °pp\ , (30) 



where I = A, B,C and 



I / U2 

e A = a , e B = a~ l , e c = -\ —p-. (3i) 

a V -H 
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Making use of the gauge condition results in the following expansions for 
F and the various components of ||, 

F = -K^z+(v B + ^jf + AF, (32) 

E 00 = V B z-K^f + ^D A h + AE°\ (33) 
E° l = ~V B v t + AE° l , (34) 



E lJ = -V A [h %J --5 %J hj+AE^. (35) 

The mixing in (|35|) between hij and its trace h = hkk can be removed by 
taking linear combinations. The same is true for the mixing in (R3|) between 



h and the variables / and z. However, the mixing between z and / in fl32 



cannot be removed algebraically for a general background. To diagonalize 
this sector we must make a nonlocal field redefinition, 

x(r],x) = z'(r],x) + -<p' (r})f(ri, x) , (36) 

When all this is done the fully diagonalized equations take the form, 

V A x = — (AE 00 + AE u y — ^(p' AF , (38) 

Vcv = -^'o (Ai? 00 + AE") - <f/ , (39) 

V BVl = 2AE°\ (40) 
V A hij = -2 (AE ij - 5ijAE kk ) . (41) 

Finally, one must check that the gauge conditions can be consistently 
imposed as constraints. Their 3 + 1 decomposition in the diagonal variables 
takes the form, 

F = -2x + v ¥ - | (^\ , Fi = ~ (avi)' - h lhj - hi ti . (42) 
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Using the "commutator" identities, 



V B ad - = ad -V A , V A -d a = -d aV B , 
a a a a 



(43) 



and the gauge-fixed field equations (|31|{il"|) , it is easy to show that the con- 
straints obey, 



V B F 
V A Fi 



2 \AE^° 



-2 I 



AE^ „ + -AE 

\ ,M a 



(44) 
(45) 



The quantities on the right hand side of (44-45) vanish as a consequence of the 
background Bianchi identities in the usual way. Therefore, the constraints 
are preserved by the gauge fixed evolution equations and full equivalence 
between the invariant field equations and those of the gauge-fixed theory will 
hold if the initial value operators are constrained to make F^ and its first 
conformal time derivative vanish at rj = i] . 



3.2 Unconstrained plane waves 

It is useful to consider plane wave solutions to the linearized, gauge-fixed field 
equations before suppressing the physical gravitons and using the constraints 
and the residual gauge freedom to purge unphysical initial value data. That 
is, we seek the kernel of the differential operators T>j — k 2 = —d^ — k 2 + 9"/9i, 
where the functions 9i{rf) were defined in relation (pl|). Although explicit 
solutions cannot be obtained for an arbitrary inflationary background, the 
limiting cases of k = and 9j = serve as the basis of series solutions for 
the "infrared" (k 2 < 0//0j) and "ultraviolet" (k 2 > 0"/0j) regimes. In the 
ultraviolet regime the particle interpretation is the same as for flat space so 
its normalization is used to define that of the infrared regime. The final step 
of this subsection is superposing plane waves multiplied by the creation and 
annihilation operators of the linearized, gauge-fixed action. The commutator 
functions between such fields give the retarded Green's functions needed to 
integrate the perturbative field equations. 

The following two solutions comprise a useful basis set in the infrared 



10 



limit (jfe = 0),Q: 

QioAv) = Oi(v) , Q2oAv) = 0i(v)£^d V '-^. (46) 

One of these is a "growing" mode, the other a "decaying" mode. For example, 
the decaying C mode is Qio t c(v)- The growing C mode can be evaluated in 
the slow roll expansion, 



(47) 



We can find the inverse of T>j on any function f(rf) by two simple inte- 
grations, 

ipj x f) (v) = -Bib) fdrf-^jr /V<WW) . (48) 

By integrating the full equation, T>iQ = k 2 Q, one obtains a relation, 

Q l Av,k) = Q l0 ,i(v) + k 2 (vj 1 Q t j)(r ] ,k) £ = 1,2, (49) 

whose iteration results in a convergent series expansion in powers oik 2 . Note 
that constancy of the Wronskian and the vanishing of corrections at rj = t]q 
allows us to evaluate it using the zeroth order solutions, 

Q2M k)Qi,i(v, k) - Qi^fo, k)Q 2 ,i(v, k) = l. (50) 

The useful basis elements for the ultraviolet limit (8j = 0) are e~ lkr> 
and its conjugate. In this limit we know the absolute normalization from 
correspondence with flat space. Coupling this with the harmonic oscillator 
Green's function gives the integrated mode equation, 

e -ikr, 1 Qltr=\ 

QAv, k ) = ^ + kJ_ 00 df i sin - ^eM Ql ^ k) ' (51) 



whose iteration yields an asymptotic expansion in powers of 1/k. The Wron- 
skian for the ultraviolet mode functions is, 

Qi(v, k)QKv, k) - QT(V, k)Qi(rj, k) = -i. (52) 



1 Here our notation differs slightly from Q in that the signs of 6c and of the second 
solutions Eq. (J4q) have been inverted. 
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Any quantum operator ipi(r],x) which is annihilated by Dj can be ex- 
pressed as a superposition of plane waves with operator coefficients, 

Mv,x) = /^{Q / (r/,fc)e*^(fc)+QK^^)e-^ t (fc)} • (53) 

If the conjugate momentum is 1/4(77,0;) then the Wronskian (|52|) implies that 
the coefficients have the algebra of canonically normalized creation and an- 
nihilation operators, 

(£), = Sij{27c) 3 5 3 {k - k') . (54) 

By taking the commutator of two such fields we obtain a sequence of lovely 
expressions for the retarded Green's function of the differential operator £>/, 

Gj(x; x') = -i9(A v ) [Mv, 5), Mv', #)] , (55) 
= -iO(Ar]) J 703^ A2 {Q/fa, *) - (77 - V)} , (56) 

—e*^{Q 1J {r,,k)Q 2 j(T/,k)-(7,~rf)} , (57) 

where Arj = Tj — rf and Ax = x — x'. The crucial final expression — in terms 
of the infrared mode functions — follows from the fact that that it obeys 
the same differential equation — T>jGi(x; x') = 5 4 (x — x') — and retarded 
boundary conditions as the first two. 

The preceding analysis suffices to define the retarded Green's functions 
needed to integrate equations ( fi~0| ) and (|4lD , 

Vi (r],x) = vf n (r],x) + 2 T drj' [ d 3 x'G B (x;x')AE 0i (x') , (58) 



hij(r),x) = h%*(r),x) 



-2 (6*5,1 - SijSkt) P drj' I d 3 x'G A (x; x')AE kl {x') . (59) 

It would be straightforward to give plane wave expansions for v\ m and h l ™. 
However, we will not bother since these linearized fields vanish when physical 
graviton degrees of freedom are suppressed and the various constraints and 
residual gauge conditions are imposed. 



12 



The same general technology can be applied to integrate the more com- 
plicated (/, z) system, 



■Wo 



V 



B 



f 

Z 



AF 



AE 00 + AE 1 



The desired retarded Green's functions, 



(60) 



-K<p 

Vb 



G ff (x;x') G fz (x;x') 



G zf (x;x') G zz (x;x' 



S 4 (x - x') 



1 
1 



(61) 



also follow from commutators of the linearized fields, 



G ff (x;x') G fz (x;x') 
G zf (x;x') G zz (x;x') 



-i9(Arj) 



lf(x),f(x>)} [f(x),z(x')] 
[z(x)J(x')) [z(x),z(x')} 



(62) 

Unfortunately it is the fields x(r), x) and y(r), x) which are annihilated by T>b 
and V c , so it is they that possess simple plane wave expansions in terms of 
the B and C mode functions. 

At linearized order the gauge-fixed action reveals the conjugate momenta 
to / and z to be simply their conformal time derivatives. Definitions (j3G|) 
and (^) give the following nonzero equal-time commutation algebra for x 
and y, 

[x(r], x),x'(r], y)\ = -iV 2 <5 3 (£ - y) = [y(rj, x),y'(r], y)\ . (63) 
If we still employ canonically normalized creation and annihilation operators, 



X(k),X(k') = (2iry5 (k - k') = Y(k),Y(k') 



(64) 



then the additional factor of —V 2 in ( |63D implies the following plane wave 
expansions, 



x(rj, x) 
y(rj,x) 



d 3 k 

d 3 k 
(2^)3 



k {Q B (v, k)e rk *X(k) + Q* B (V, k)e' lts X\k)} , (65) 
k{Q c {v,k)e its Y(k)+Q* c (r 1 ,k)e- it£ Y\k)} . (66) 
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Since the transformation between (/, z) and (x, y) involves conformal time 
derivatives, its inverse cannot be local in time for the off shell fields. However, 
for the on-shell solutions we can use the linearized equations to obtain the 
following expressions, 



/ = 



z = 



V 2 
1 



-ip' x + y' + ^y 



x - 



(67) 
(68) 



It is simplest to exploit these relations in evaluating the various commutators 
of the (/, z) retarded Greens functions. One then replaces products of the 
mode functions with the B and C-type Green's functions. The answer is, 

„2 



G ff (x;x') = ^ <j -— w' {ri)w' Q {ri')GB{x;x') 



G fz (x;x') = 
G zf (x;x') = 



K 



- do&MvWMGcix; x>) + 5\x - x') 
Vo(v)9' Gb{x; x') + ^do<p' (v)Gc(x-,x')} 
<p' (rj')d G B (x;x') + ^f'M)G c (x- x>) 



2V 2 

K 

2V 2 



(69) 
(70) 
(71) 



G Z z {x, X ) 

1 

V 2 



-8 & Q G B {x- x') - ^ V ' ( V )v' W)G c (x; x') + 5\x - x') \ .(72) 



Finally, we substitute the relation, Gi(x; x') = V T 5 4 (x — x'), 



Gf f (x;x') 
G fz (x;x') 
G zf (x;x') 

G zz (X) x ) 



1 

V 2 

K 

2V 2 

K 

2V 2 
1 

V 2 



K 2 1 "! 

4 ^0 



5\x - x') . (73) 



-(p'oV^do + — doip'oV^ip'o 



1 + doV^do - ^Vo^cVo 



5\x - x') , 
5\x - x') , 
5\x - x 1 ) , 



(74) 
(75) 
(76) 
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3.3 Integrated field equations 

In this subsection we first purge unphysical initial value data by imposing the 
constraints and making use of the residual gauge freedom. Then gravitons 
are suppressed to leave what we call the physical linearized fields. Finally, the 
various field equations are integrated using the retarded Green's functions of 
the previous section. 

Since our gauge condition involves derivatives it is preserved by residual 
gauge transformations which obey second order differential equations. 2 One 
can therefore freely impose four residual conditions and their first conformal 
time derivatives on the initial value surface. We choose these to null x(r] , x), 
x'(r]o,x), Vi(r} ,x) and v'^^x). 

At the end of subsection 3.1 we showed that the four constraints are 
also annihilated by second order differential operators. It therefore follows 
that the initial values of the four constraints and their first conformal time 
derivatives are the only extra conditions to be imposed in addition to the 
gauge-fixed field equations. With our residual gauge choice these eight con- 
ditions imply that hij(r),x) is transverse and traceless at linearized order. 
Suppressing gravitons makes it vanish completely at linearized order. 

The only linearized fields which remain are / and z. Since x vanishes, 
they can depend only upon the operators Y(k) and Y*(k), 

/ P hM = _^ (^)y_-{g o (^A ; ) c <,Mf y(^)+h.c.} . (77) 

ZphM = ^o(v) J ^l{Qc(v,k)e rk x Y(k) + h.c} , (78) 

We can obviously eliminate / p h in favor of z p ^, 

f , - 2z' h (r],x) z' h {r],x) 

fphiv, x) = = — r = . 79 

This is the well-known constraint between the (dimensionless) Newtonian 
potential $(77, x) = Kz(rj,x)/2a(rj) and the scalar field fluctuation [11 . 

We can now write down the integrated field equations whose iteration 
produces the final perturbative expansion of the various fields, 

f(r],x) = /phfaa!)- fdrf [ d 3 x' \G ff (x;x')AF{x') 



The analysis is quite similar to that carried out for de Sitter background in [[l0| 
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+ G fz (x; x') \AE 00 (x') + AS"(a/)l \ , ( 80 ) 



z{rj,x) 



z ph (r],x) - [ drf [ Sx \G zf (x; x')AF(x') 




Vi(rj,x) 



2 f drf I d 3 x'G B (x;x')AE 0i (x') , (82) 



hij(rj,x) 



-2 (5 lk 5 je - 5 {j 5 M ) [" drl [ d 3 x'G A {x; x')AE k \x') . (83) 



Note that Vi and h%j vanish at rj = tjq, along with their first derivatives. On 
the other hand, both / and z suffer perturbative correction on the initial 
value surface. This derives from the existence of a gravitational interaction 
even on the initial value surface, as required by the constraint equations of 
the ungauged formalism. 

3.4 The long wavelength approximation 

Except for suppressing the gravitons, all the results obtained to this point 
have been exact. Unfortunately they are also largely useless because we lack 
explicit forms for the various mode functions and retarded Green's functions 
in a general inflationary background. Nor could we perform the required 
integrations if we did possess such expressions. What makes calculations 
possible is the long wavelength approximation which effectively removes spa- 
tial variation. The resulting problem of temporal calculus is still formidable, 
but it can be treated using the slow roll expansion discussed in Section 2. 

The back-reaction we seek to study is the response to the superadiabatic 
amplification of modes which experience horizon crossing. This is an in- 
trinsicly infrared phenomenon so one might suspect that it can be studied 
effectively in the limit where spatial derivatives are dropped and the mode 
functions approach their leading infrared forms. This limit defines the long 
wavelength approximation and it effects three sorts of simplifications on the 
perturbative apparatus of relations ( ^0 - ^3 ) . 

The first simplification is to ignore spatial derivatives in the source terms, 
AF and AE^ V . An immediate consequence of this, and of the absence of 
dynamical gravitons, is that AE 0t vanishes and AE %3 is proportional its 
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trace, 



A£«|,„.^0 , A^W^-A^. (84) 

It follows that Vi remains zero for all time and that only the trace part of hij 
ever becomes nonzero, 

Vi(r],x)\ e . w . — > 0, (85) 
h l3 (ri,x)\ Lw . — > % fdrf f d 3 x'G A (x;x')AE kk (x') . (86) 

The second simplification is that one can replace the mode functions of 
the linearized solutions by their infrared limiting forms. These forms were 
worked out in J3J by matching the leading ultraviolet term of the ultraviolet 
expansion to the leading term of the infrared expansion at the time r]*(k) of 
horizon crossing: k = H(r)*)a(r)*). Since the physical fields contain only C 
modes the result we require is, 



Qc(v, k)U.w. — ► \ 2( v 1 . —,=■ ■ 87 



Note also that momentum integrations are cut off at k — H(7])a(i]). The 
physical justification for this is that only modes below the cutoff have under- 
gone the superadiabatic amplification that is the basis for the effect whose 
back-reaction we seek to evaluate. However, one must take care that the re- 
sulting time dependence is a genuine infrared effect and not simply the result 
of allowing more and more modes to contribute to what would be an ultravi- 
olet divergence without the cutoff. One consequence is that time derivatives 
of the fields should not be allowed to act on the momentum cutoff. 

The third simplification is that we expand the inverses of T>i = T>j + V 2 
in powers of V 2 , 

VJ 1 = VJ 1 - VJ 2 V 2 + ... (88) 
Only the zeroth order term is required for h^, 

h i3 \ Lw . — > ^6 ij Vl 1 AE kk , (89) 

but one must go to first order for the various (/, z) Green's functions on 
account of their prefactors of 1/V 2 . The Appendix demonstrates that the 
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leading results are, 



Gf f (x; x')\ Lw . 
Gf z (x; x )\tw. 
G zf (x;x')\ tw . 



G zz {xy x ) 



V^o </o V^o 

k^V^^VcVo^x-x'), 

k^Vc^VI V 5 4 (x - x') , 

(f ' V c 1 \ + 2^V c 1 ^V B 1 d 



<5 4 (x-a/) 



5\x-x') ,(90) 
(91) 
(92) 
(93) 



At this point we can make contact with the similar calculation that was 
done in In the infrared limit, the general retarded propagators (j90"|)-(j93|) 
correspond exactly to the result that can be inferred from Eqs. (131)-(132) 
of J! . The "C" and U D" in Eqs. (131) and (132) from that paper (where 
Q = a) correspond respectively to the change induced in the z and in the / 
fields that obtains from our equations (|80| - |81|) . The the amputated 1-point 
functions in those equations correspond to our source terms according to the 
rule: a -> -f aAE k \ 7 -> —naAE 00 and 6 -> -naAF. 



4 Quadratic corrections 

The purpose of this section is to apply the technology of Section 3 to obtain 
the scalar and metric at quadratic order in the initial value operators Y{k) 
and Y'(k), and to leading order in the slow roll and long wavelength expan- 
sions. In Ref. we computed the one loop expectation values of the scalar 
and metric fields, also to leading order in the slow roll and long wavelength 
expansions. Much of that calculation can be used to get the quadratic order 
operator expansions of the same quantities. In particular, all the vertices are 
catalogued in that reference. 

As an example, consider vertex # 9 of Table 3 from |4]]. That vertex 
corresponds to the following cubic term in the Lagrangian, 

£#9 = aW = -~V tVV (<Po) + h) . (94) 
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Variation with respect to 4> an d gives the following source terms 



AF #9 = -KV tVV (<po) af(z + -h) 



^(^)a/V v . (95) 



The lowest order perturbative corrections from these terms come when /, z 
and h are replaced by the associated physical linearized fields: / p h, z p h and 

0, respectively. Recall also that / ph = z'/ayJ—H. 

The preceding are all exact results. Because we are also making the long 
wavelength and slow roll approximations, z p h can be written as, 



z ph (r],x)\t 



w.+s.r. 



a{r)) 



d 3 k H* 



(2tt) 3 -H, 



m 



* F {Y(k) + Y^(k)} 



(96) 

To leading order in the slow roll approximation the time depedence of this ex- 
pression comes entirely in the initial factor of a(rj), which allows the following 
simplification for 



w.+s.r 



z ph (r},x) . 



(97) 



We can also use the slow roll approximation to evaluate derivatives of the 
scalar potential, 



6H\ -H 



Vipip(<fo)\s.r. 



-3H . 



(98) 



The source terms finally contributed by vertex ^9 in the slow roll and long 
wavelength approximations are, 



AF 



#9\Lw.+s.r. 



A /? 00 I 

'-*- c '#g\Lw.+8.r. 
A TT'kk I 

'-± I1/ #9\£.w.+s.r. 



-H 



JJ2 aZ ph ' 



+ -KH 2 azl h , 
g 

— ~kH az ph . 



(99) 
(100) 
(101) 



Summing the contributions from all vertices of Einstein-scalar Q under the 
same set of approximations gives the following quadratic source terms, 



AFI 



l.w.+s.r. 



3kH 2 



-H 



fj2 aZ ph ' 



(102) 
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AE 00 \ e . w . +s . r . — ■* -~KH 2 az* h , (103) 

27 

A£ fe \„, +S . r . — +— wtf 2 ^. (104) 

The next step is to apply the various retarded Green's functions in the 
long wavelength and slow roll approximations. As an example, we substitute 
( |104|) into (|89|) to obtain the leading quadratic correction to hy, 

h-ij \tw.+s.r. — > — 9/c5y-a / dr/' ' oT 2 I dr/" a 2 H 2 z 2 h , (105) 

Jt]o J "no 

where z p h is given by Eq. (|96|) . Taking the momentum integrals outside, we 
are left with the following integrations over time: 

a(t) £ dt'a-\t>) £' dt"a\t")^- 2 = ~a(t) (^j +... (106) 



The result is therefore, 

/ H 2 \ f 3 kz 2 



h l3 \ Lw . +s . r . — ( ^ W + 0(z^ h ) . (107) 



The analogous reductions for / and z are straightforward but rather tedious, 
owing to their complicated mixing. The final answer is, 



Up { 5 m 2 



f\t w .+s,. — fjp + 4V+%)j' ( 108 ) 
2k«,.+.. r . — ► {^ph + ^^ + C(4h)} • (109) 



We invert (|29|) to recover and ^j,. It is useful to multiply by k to absorb 
the dimensions, and to express the result in terms of the (dimensionless) 
Newtonian potential: 

^ KZ^ 

Note that, from Eqs. ( fLlOj) and fl9"6|) we have $ cx i/~ 2 (since H ~ constant) 
so that its time derivative is down from if $ by a slow roll parameter, 

$k w .+s.r. — 2 ( ^ J tf$ . (Ill) 
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With this terminology, the nonzero perturbed fields are, 



/ H 2 

K<f>\tw.+..r. — > W-gr{2$ + 5<l> 2 + 0($ 3 )} . (112) 

«^ook™. s .r. — > 2$ + 14$ 2 + 0($ 3 ) , (113) 

«^y|/.«,.+..r. — > ^ |2$ - 6 f-^T) $ 2 + 0($ 3 )| . (114) 

5 An invariant measure of expansion 

The purpose of this section is to compute corrections through quadratic order 
to the operator A im/ we have proposed as an invariant measure of the rate 
of cosmological expansion. We first expand the scalar A = D^ 1 in powers of 



the pseudo-graviton field and then substitute (|113[ - |114D to express the result 
as a function of the unconstrained initial value operators, to leading order in 
the slow roll and long wavelength approximations. Full invariance is achieved 
by evaluating the A(r), x) in a geometrically specified coordinate system. 

5.1 The scalar observable 

The pseudo-graviton expansion is most easily accomplished by first express- 
ing D c in terms of the conformally rescaled metric, 

g^(v, x) = a~ 2 {ri)g^{ri } x) = rj^ + K^fa x) . (115) 

We write O c = a~ 3 Da, where T> and its expansion are, 

1 „ / r-i_..,.„ \ 1 



V = [\f-9 <TUj -~R = 2 + kPx + k 2 V 2 + ... . (116) 

The first two operators in the expansion are, 

v 1 = -r v d»d v + (-r« a + \^) p v - > ( 117 ) 
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The second order, conformally rescaled Ricci scalar is, 



4 f~> 2 

The next step is to factor d 2 out of V, 

V = d 2 (l + ^kD x + ^k 2 V 2 + 0(k 3 )) . (120) 

Inverting T> is now straightforward, 

111 11 1 1 1 , 1 „ , , 

V = ¥ ~ V^W 2 + & KVl & KVl & ~ ~~9^ 'd 2 + ° {K ) ■ (121) 

All this implies the following expansion for the scalar observable, 

A[g] = a-^a 3 = A + kAi + k 2 A 2 + 0(k 3 ) , (122) 
where the first two corrections are, 

A 1 = -ar^V^a\ (123) 
= -^^V^ + a^^V^V^a*. (124) 

The zeroth-order term Aq = a d~ 2 a s has a very simple expression in 
terms of the background expansion rate, 

* = £ jC * W) = 2^ { 1 + ° (i^) } • (125) 

Our proposal is to define the general expansion rate, as an operator, to bear 
the same relation to the full scalar A |J. Astronomers do the same thing 
when they infer the Hubble constant under the assumption that the relation 
between luminosity distance and redshift is the same in the actual universe 
as for a perfectly homogeneous and isotropic one. 

The higher order contributions A± and A 2 are nominally inhomogeneous 
but lose their dependence upon space in the long wavelength approximation. 
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The consequent suppression of spatial derivatives and anisotropic compo- 
nents of the pseudo-graviton field reduces T>\ to the following simple form, 

V x \ Lw , — > -^oo9 2 - + iQdo - ~V4 . (126) 

Further simplifications result from the fact that (co-moving) time derivatives 
of $ are weaker than if $ by a slow roll parameter as per equation fllllj) . Only 
the quadratic term of ijjy possesses the enhancement necessary to survive 
differentiation at leading order, 

^Utw.+s,. — ► #a{o + 0($ 3 )} , (127) 
Kiff'ijk^.+s.r. — > #a{-12$ 2 + 0($ 3 )}c% . (128) 

Of course subsequent conformal time derivatives are dominated by the fac- 
tor(s) of a, 

K1pwU. W .+s.r. — > #V{0 + O($ 3 )} , (129) 

k^Ilw.+s.t. — > #V{-12$ 2 + 0($ 3 )}c% . (130) 

The result is that only a few terms in Pi and T> 2 can contribute at leading 
order, 

«£>ik™.+ s .r. — ► -2$3 2 + $ 2 {-14<9 2 + 18iia«% + 6F 2 a 2 } , (131) 
«P 2 k. w .+ s .r. — > -4$ 2 <9 2 + 0($ 3 ) . (132) 

Because $ <C ff$ we can ignore its time dependence with respect to 
factors of a(i]). This makes acting the various factors of d~ 2 quite simple in 
the slow roll approximation, 

Axku.+s.r. — > -tU 3 (-2$ -2$ 2 + 0($ 3 )) , (133) 



a <9 2 

— > ^{-2$-2$ 2 + 0($ 3 )} . (134) 
Since the two terms in A2 cancel to leading order, the result is, 

A(x)\ Lw . +s . r . — > ^ {l - 2$ - 2$ 2 + 0($ 3 )} . (135) 
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5.2 The invariant observable 



Even a scalar changes under coordinate transformation. To achieve full in- 
variance we measure ^4.[g](x) at a point Y^[(f, g](x) which is invariantly speci- 
fied in terms of the operators ip and g^ v . Since there is no spatial dependence 
in the long wavelength approximation we actually need only fix the 3-surface, 
Y° = t. Following the suggestion of Unruh, we chose the functional r[(p] (rj, x) 
to make the full scalar agree with its background value at conformal time rj, 



<p(r{r],x),x \ = <p (ri) . 



(136) 



The weak field expansion of r[ip] is 



r[ip](r},x) = r}- 



(j)(r],x) 6(rj,x)(j)'(rj,x) (p'^(rj) (<f>{r), x ) 



<Po(v) 



2<Po(v) V <Po(v) 



+O(0 3 ) . (137) 



Substitution of our result ( |112| ) for the perturbed scalar gives the final expan- 
sion to leading order in the slow roll and long wavelength approximations, 



7] + 



H 2 



-H / Ha 



H 



2 \ 2 



7 



H 2 



2Ha 



+ 0($ 3 ) . (138) 



-H) \-H / 
Our invariant expansion operator is just A evaluated at this point, 

Anvb, #](??,£) = A[g) { r[<f)(r],x),x ) . (139) 



Note that there is no obstacle in perturbation theory to evaluating an oper- 
ator at a point r[(p] = rj + 6r[ip] which is itself an operator, 



A mv = A+ A'5t + 2>A"5t 2 + ... 



The derivatives are straightforward to evaluate with (|135|) and ( |1 1 1|) , 

-Ha ( -H' 



A \l.w. J r s.r. 

A"\ 

<s\ \£.w.+s.r. 



2H 2 \H 2 
-H 2 a 2 



{2-8$ + 0($ 2 )} 



-H 



-H 



2H 2 



H 2 



H 2 



+ 0($) 



(140) 

(141) 
(142) 



24 



However, putting everything together results in complete cancellation of all 
corrections to the order we are working, 



A„v = ^{l + 0($ 3 )} • (143) 

Although this result surprised us, it could have been anticipated by not- 
ing that, in the slow roll approximation, our scalar degenerates to a local 
algebraic function of the Ricci scalar || , 

AWWU^j^- ( 144 ) 

Einstein's equations completely determine the Ricci scalar, as an operator, 
from the local matter stress tensor, 

R(x) = -SnGg^T^, (145) 
= IGttG - X -g^d^d v ^ . (146) 

We can always choose to work in a coordinate system for which the full 
scalar agrees with its background value. When this is done one sees that 
back-reaction can only enter through the kinetic term. Since it is certainly 



from kinetic effects that the pure gravitational result derives |[12j| , there may 
well be a significant back-reaction from the scalar kinetic term as well. Un- 
fortunately, one can never tell whether there is or not when the kinetic term 
is systematically neglected, which is just what happens when the slow roll 
and long wavelength approximations are combined. It follows that our ap- 
proximations must produce a null result — whether or not there really is 
significant back-reaction — not just at quadratic order in the initial value 
operators, but at all higher orders as well. 



6 Stochastic samples 

The last section has demonstrated that scalar-driven inflation can show no 
secular back-reaction to leading order in the long wavelength and slow roll 
approximations. Since this holds as a strong operator equation for our ex- 
pansion observable, A im/) there is no need to choose between expectation 
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values and stochastic samples. However, we believe the case is still quite 
strong for secular back-reaction when the long wavelength approximation is 
relaxed, although it would have to come from the coherent superposition 
of interactions at higher than quadratic order in the initial value operators. 
The purpose of this section is to develop a theoretical framework for study- 
ing such an effect through stochastic samples as recommended by Linde and 
others f7]. We begin by motivating and defining the basic formalism, then 
we treat the crucial question of the degree of stochastic fluctuation expected 
in functionals of the fields such as Anv[^, <?]• 

6.1 Motivation and basic formalism 

The simplest way to motivate stochastic effects is by considering the linearly 
independent mode functions of the infrared regime, Qij(r], k) and Q2,iij]i 
which were described in subsection 3.2. The full linearized field operator 
must involve both of these mode functions in order to avoid commuting 
with its conjugate momentum. However, after horizon crossing (k = Ha) 
one of the mode functions becomes vastly larger than the other. This is 
the phenomenon of superadiabatic amplification and it corresponds to the 
simple physical picture of particle production through infrared virtual quanta 
becoming trapped in the inflationary Hubble flow. If one retains just the 
larger mode function then the linearized operator which results is effectively 
classical in that it commutes with its time derivative. Although the result 
is still probabilistic, one can simultaneously measure the value of such an 
operator and its time derivative, and subsequent measurements will show 
only the classical evolution expected from these initial values. This is the 
reflection, in the Heisenberg picture, of a "squeezed state." 

The expectation value of a functional of squeezed operators can fail to 
provide a good estimate of what an actual observer sees. For example, the 
expectation value of the current stress tensor is presumably homogeneous 
and isotropic. We see nothing of the sort because we exist at the end of 
a long period of essentially classical evolution from one particular choice, 
random but definite, of superadiabatically amplified density perturbations. 

A better way of treating squeezed operators is to sample the result of 
making random but definite choices for them from the relevant quantum 
mechanical wave function, and then evolving classically. It is crucial to un- 
derstand that taking such a stochastic sample is perfectly consistent with 
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the use of quantum field theory to express the Heisenberg operators as func- 
tionals of the unconstrained initial value operators. Nor is there any change 
in how the observable Ai nv (r], x) depends upon the Heisenberg operators. 
(Of course we do want to avoid making the long wavelength approximation!) 
What changes is that the scalar creation and annihilation operators — Y\k) 
and Y(k) — are random numbers up to k = H(j])a(j]), and zero beyond. 

To avoid problems with continuum normalization we take the 3-manifold 
to be T 3 with identical co-moving coordinate radii of Hq 1 = H^ 1 ^). (Be- 
cause the conformal coordinate volume is so restricted during inflation the 
integral approximation is excellent for mode sums and there is no conflict 
with any of our previous, continuously normalized, results.) On this mani- 
fold the co-moving wavevectors become discrete, 

k = 2nH n . (147) 

Phase space integrals are converted into mode sums in the usual way, 

d 3 k 



/(I fc -* 
J ^ V3 f(k)^H^f(2nH n). (148) 



{2<Kf 

And the Dirac delta function goes into a Kronecker one, 

(2tt) 3 5 3 (£ - k') — -> H 3 5 HiH , . (149) 



Because the initial scalar state is free, the associated creation and anni- 
hilation operators are stochastically realized as independent, complex, Gaus- 



sian random variables with standard deviation H n 3 . It is convenient to scale 



out the dimensions, 

A H = HlY(2nH n) , ^ = #j ^(27^) , (150) 
so that the probability density for each mode n has a simple expression, 

p(An,A* n ) = ±e- A * A *. (151) 

Since the modes are independent, the joint probability distribution is just 
the product of (|151|) over the relevant range of n. 
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6.2 Functionals of stochastic variables 



Our observable Ai nv (r], x) depends in a complicated way upon the scalar and 
metric fields, which are themselves functionals of the stochastic variables 
Afi and At. The nature of this dependence determines the crucial issues of 
whether or not Am V (rj, x) is well represented by its expectation value and 
whether or not a definite sign can be inferred for corrections to the cosmo- 
logical expansion rate. Of course we do not yet have a replacement for the 
long wavelength approximation which shows secular back-reaction. However, 
we do here develop a technique for separating nonlinear functionals of the 
stochastic variables into a part whose percentage fluctuation becomes negli- 
gible for a long period of inflation, plus another part whose fluctuation is not 
negligible but which has a definite sign. 

Recall from probability theory that a functional f = F[A, A*} of random 
numbers is itself a random number. Its probability distribution function 
descends from that of the A^ and At by Fourier transformation, 

p(f) = JZ^ eikf ( e ~ ikF[A ' A * ] ) ' (i52) 

= r ^ e */ (jj f M*^k e -A*A$\ e -ikF[A,A*] . (153) 
J—oo 2tx y ^ i 2lX J 

As an example, consider the Newtonian potential ( |HU| ), 



*(JI, *) = ^ H " H{r,) E - {Qdv, 2*H Q n)^ H ^A n + c.c.} . (154) 



ri 



The various integrations are trivial Gaussians, as they would be for any linear 
variable. The result is that $ follows a Gaussian distribution with mean zero, 

m = -^= 2 ^i^ , (155) 

and a spatially constant variance equal to, 

° 2 (V) = Z ^ L E^\\Qc(v,2nH n)r . (156) 

ri 
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Making the slow roll and long wavelength approximations we see that the 
variance grows as the logarithm of the scale factor, 



The simplest sort of nonlinear functional is just the square of a linear 
one. Since linear functionals of A^ and At are always Gaussian, their squares 
always follow a Chi-squared distribution whose mean is the variance of the 
Gaussian and whose variance is twice the square of this. For example, the 
variable $ 2 (^, x) follows a Chi-squared with mean a 2 and standard deviation 
\f2o 2 . Although the fluctuations of $ 2 are of the same order as its mean, 
the sign is definite. Note also that, for a long period of inflation, only an 
incredibly fortuitous sequence of choices for the stochastic variables A^ and 
At would result in $ 2 having a value significantly below some constant times 
our estimate ( |159|) . So that if some reliable approximation scheme should 
wind up giving the effective expansion rate as, 



then the conclusion would be that there is only a vanishingly small proba- 
bility to observe anything except a secular slowing of inflation. This is an 
example of how stochastic samples might show the same qualitative results 
as expectation values while breaking exact homogeneity and isotropy and 
altering the numerical coefficient of the order of magnitude estimate. 

It might be thought that quantitative results are not obtainable for more 
complicated nonlinear functionals. That useful statements can still be made 
derives from the following fact: any secular back-reaction effect must result 
from the coherent superposition of contributions from an enormously large 
number of modes. This allows us to exploit the same sorts of simplifications 
that underlie statistical mechanics. 

To illustrate the important considerations without becoming too mired 



o- (T])\lw.+ S .r. 




H eS (r],x) = H(r]) {l + $(77, f) - $ 2 (7?,£)} 



(160) 
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in technical detail let us consider quadratic superpositions of the form, 



s — S[A, A ] = — 2J ^mn (Art + A^) (Afi + A H ) 



The characteristic function is, 

-ikS[A,A*] 



det{I + 2ikS) 



(161) 



(162) 



where S stands for the symmetric matrix Sjnn and I is the unit matrix of the 
same rank. The various moments of s follow by differentiation, 



(s) = Tr[S] , i(s 
Tt[S}) 3 ) = 4Tr[S 3 ] , 



Tr[S]) 2 ) = 2Tr[S 2 



(163) 
(164) 

How much fluctuation one should expect is governed by the relation be- 
tween traces of powers of the matrix S. We distinguish two cases: 

1. "Local" superpositions which are characterized by Tr[S ,?1 ] = (TrfS 1 ])" 
for n > 2; and 

2. "Nonlocal" superpositions which obey TrfS 2 ] ^ (Ti[S]) 2 and TrfS" 1 ] 
Tr[S 2 ](Tr[,S]) n ~ 2 for n > 3. 

The local case is just the square of a linear superposition and has already 
been considered. To recapitulate, it has significant fluctuation but definite 
sign. The distribution for a nonlocal superposition can be approximated by 
dropping higher traces of S in the exponent of ( |153] ), 

f°° dk 



Pis) 



= ifcs-iTr[ln(/+2ifcS)] 



-00 27T 

00 dk 
— ( 
-00 2n 

1 



t iks-ikT[[S]-k 2 TtlS 2 



AttTt[S 2 } 



exp 



Tr[S)y 



4Tr[5 2 



(165) 
(166) 
(167) 



This is just a Gaussian centered on Tr[S] with variance 2Tr[S' 2 ]. Since its 
standard deviation is insignificant compared with the mean, stochastic effects 
are not important. 
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Recall that invariant measures of the cosmological expansion rate derive 
their dependence upon and A* H from two expansions. In the first the ob- 
servable is expanded in powers of the perturbed fields, then one substitutes 
the expansions of the perturbed fields in powers of the stochastic initial value 
data. Both of these expansions really involve averages over the past lightcone 
of the observation point. However, depending upon the observable and the 
approximation techniques used to solve for the fields, the terrific expansion 
of spacetime during inflation may weight the averages heavily towards the 
most recent times. In that case one gets a local superposition and stochastic 
effects are important but simple to include. An improved observable can also 
be defined so that the entire past lightcone participates effectively 0, and it 
seems likely that improving on the long wavelength approximation will have 
this result anyway. In this case one tends to get a nonlocal superposition 
because modes with different n interfere destructively. The stochastic fluc- 
tuation of back-reaction would then be negligible and one may as well resort 
to expectation values. 

There is also the possibility of mixing of local and nonlocal, in which 
case the best strategy is to separate the two effects and treat them as above. 
One might anticipate that analytic considerations would render the local part 
obvious but even if not, its dyadic form makes the decomposition a simple 
linear algebra problem. For example, suppose that at late times we have, 

Tr[5 n ] — ► («Tr[5]) n , (168) 

for some positive constant a < 1. The putative decomposition would be, 

Sfhn "mUn ~\~ ASVjm j (169) 

where the various traces obey, 

u-u — ► aTi[S] , Tr[AS] — ► (1 - a)Ti\S] , Tr[AS n ] < (Tr[S]) n . 

(170) 

Simply pick any mode function with nonzero overlap — it need not be 
close to tin — then contract into S 2 and divide by the trace of S, 

(S 2 v)h u*ASv u-v, An . (AS 2 v) H 

Wr = au ' " "* + "w " fl + ^si {ASuh + ^jsr • (171) 

Only the first term can matter at late times so we recover ua by normalizing 
and then multiplying by the square root of aTrfS']. 
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7 Summary and discussion 



We have calculated the gravitational back reaction on scalar-driven infla- 
tion using an invariant observable, to quadratic order in the initial creation 
and annihilation operators and to leading order in the long wavelength and 
slow roll approximations. No effect was found, contrary to previous work 
by ourselves and others which indicated a secular slowing of the expansion 
rate at this order [§, [|, f|]. It is significant that the inclusion of stochas- 
tic effects played no role in this change. Under our approximations, secular 
back-reaction would enter through the accumulated Newtonian poten- 

tial from modes which have redshifted into the infrared regime. However, all 
terms of order $ and <3> 2 drop out of the Heisenberg operator Ai nv , before 
one has to choose between the alternatives of expectation values or stochastic 
samples. 

Another improvement in our analysis which had nothing to do with chang- 
ing the result was the use of a scalar to measure the expansion rate. The 
entire difference from previous work is in fact attributable to defining the 
coordinate system so as to make the quantum scalar vanish on surfaces of 
simultaneity. When these coordinates are used even the gauge fixed expec- 
tation value of the metric fails to show significant back-reaction at one loop. 

The null result was anticipated by Unruh who noticed that scalar mode 
solutions become pure gauge in the long wavelength limit ||. This does not 
preclude back-reaction but it does rule out dependence upon $ which fails to 
vanish in the long wavelength limit. For example, spatial and certain tem- 
poral derivatives of can contribute. In fact there are such contributions 
but they are negligible at quadratic order. 

Although we have obtained a discouraging result about the possibility of 
a simple one-loop effect, our analysis does not invalidate the idea of gravita- 
tional back reaction on inflation. It would be fairer to say that what we have 
learned constrains the form any such effect can take. In particular one should 
not take the long wavelength limit. It is highly significant, in this regard, 
that the purely gravitational effect claimed at two loops |12| does not involve 
the long wavelength approximation in any way because the locally de Sitter 
background is simple enough that the full propagator can be worked out. 
In fact the gravitational response to the inflationary production of gravitons 
comes entirely from the graviton kinetic energy and would vanish in the long 
wavelength limit. 
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The long wavelength approximation was also avoided in the effect claimed 
at three loops for massless, minimally coupled ip 4 theory in a locally de Sitter 
background []TB| . In a subsequent paper |TJ| we resolve the issue of coincident 
propagators in this model by a procedure of covariant normal ordering. The 
resulting theory exhibits a secular back-reaction which slows inflation in a 
manner that is unaltered by either the use of an invariant operator to measure 
expansion or by the inclusion of stochastic effects. 

Although there was no need to consider stochastic effects in the present 
work, Section 6 describes a formalism for including them in higher order 
processes which may show secular back-reaction. This formulation differs 
from the standard one || in three ways. First, the focus is perturbative and 
local whereas previous previous treatments have been concerned with non- 
perturbative effects on the global geometry. A second difference is that past 
treatments incorporated stochastic degrees of freedom only as they experi- 
enced superadiabatic amplification. Although this is doubtless an excellent 
approximation for the sorts of global and nonperturbative issues that were be- 
ing studied, we cannot afford to ignore the non-conservation of stress-energy 
implicit in continually injecting new degrees of freedom into the system. 
Therefore, the creation and annihilation operators for any mode we wish to 
treat stochastically are considered to be nonzero random numbers even on 
the initial value surface. The final difference is that we enforce all the pertur- 
bative field equations of the Einstein-scalar system so that the various modes 
are in gravitational interaction even on the initial value surface. 

An amusing consequence of all this is that stochastic effects provide a 
nonperturbative proof that back-reaction must eventually become significant 
if nothing else stops inflation first. For it will be noted that, in our ver- 
sion, including stochastic effects at some time rj corresponds to a universe 
that began inflation at rjo with a random collection of modes excited up to 
co-moving wave number k = H{rf)a{rj). The hypothesis that inflation never 
slows amounts to the assumption that inflation can begin with the inital state 
populated to arbitrarily high wave number. Of course this is nonsense. Even 
if one subtracts off the spatially averaged energy density — as should prob- 
ably be done — what must actually happen for arbitrarily high excitations 
is that random inhomogeneities produce a gravitational collapse. 

Finally, we comment on the degree to which a random stochastic sample 
should be expected to differ from its mean. Because secular back-reaction 
must manifest itself through the coherent superposition of an enormously 
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large number of independent random variables, one can sometimes employ 
the methods of statistical mechanics. The exception is when back-reaction 
involves an ordinary function of the local stochastic fields, for example Q 2 (x). 
This can happen if the terrific inflationary expansion causes the average over 
the past lightcone to be weighted so that only interactions just before the 
observation contribute effectively. In that case stochastic fluctuation is not 
negligible, but the sign of the effect is definite. In the other, "nonlocal" case, 
stochastic fluctuation is negligible compared with the mean effect and one 
may as well use expectation values. Our suspicion is that improving upon the 
long wavelength approximation and improving on the expansion observable 
will result in secular back-reaction of the nonlocal sort. However, if there 
should be local mixing, it is straightforward to untangle. 
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Appendix 

In this Appendix we derive the infrared limits of the generalized Green's 
functions of Section 3.4. First notice that from Eq. ( |5DD and the background 
identities we can derive the following useful "commutation" relations: 

d V B -V B d = (172) 
(p , V B -V C ip' = 2ip' d < 4 , <p' V c - VB<p'o = 2^-do<p' . (173) 
l r V B -V c \ = -2^d , \v c -V B \ = -2d ^2 . (174) 



34 



Using Eqs. ( |30"D and ([H]) for T>b and T>q we have also the relations 

V B -^do<p' = d \v c -2^V 2 , (175) 

V c \d = <p' d -^V B + 2^V 2 . (176) 

Dropping spatial derivatives and inverting implies the following "inverse com- 
mutation" relations, 



T d <f/ V o 1 (f/ = tp'oV^do 



\d V B l v ' = VaV d \ , (177) 



Substituting ( 177|) into the right-hand-side of Eq. (ff5[) we see that the 
1/V 2 term cancels because, 



K 2 Ik 2 
! (p' 4 



^ K 2 , 2 

-P C + —if/* 



(179) 



Similar manipulations reveal that the 1/V 2 terms cancel as well in the ex- 
pansions (ff3|)-(|75|). 

The next order terms can be found by expanding Vj 1 in powers of the 
Laplacian. Doing this in expression (|76|) for G zz (x;x r ) yields, 



(~*zz(3'; X ) 



-doV^V^do + - v ' V c l V c l ip' Q + 0(V 2 ) 



6\x - x') . (180) 



The term of order V° can be simplified with the commutation relations (|172 
T7§, 



-d V B L V B L d + -<p' V c L V c L <p' 



-^Vo^o— + -r^c 1 — 7r¥>o ^b^o , (181) 

(f/ Q 4 (^o »o 

-^P-V^o^^^o + V^ — iVB + SftWdo , (182) 
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+ trfPc^doV^do , (183) 

= ^Vc 1 ^ + 2^ 1 ^- 2 V B l d , (184) 

The last line, acting of a delta function, gives the long wavelength limit of the 
zz retarded propagator, Eq. (|9"3"1). Similar reductions pertain for the other 
Green's functions of the (f,z) sector, (]90"l)- (j92"]) . 
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